KIAS-P13014 
IPPP/13/14 
DCPT/13/28 



Multi-component Dark Matters in Two Loop Neutrino Model 

Yuji Kajiyama, a Hiroshi Okada, fe Takashi Toma c lf| 

a Akita Highschool, Tegata-N akadai 1, Akita, 010-0851, Japan 
b School of Physics, KIAS, Seoul 130-722, Korea 
'Institute for Particle Physics Phenomenology University of Durham, Durham DH1 3LE, UK 



Abstract 

We construct a loop induced seesaw model in a TeV scale theory with gauged U(1)b-l 
symmetry. Light neutrino masses are generated at two-loop level and right-handed neutrinos 
also obtain their masses by one-loop effect. Multi-component Dark Matters (DMs) are included 
in our model due to the remnant discrete symmetry after the B — L symmetry breaking and 
the Z2 parity which is originally imposed to the model. We investigate the multi-component 
DM properties, in which we have two fermionic DMs with different mass scales, O(10) GeV 
and O (100-1000) GeV. The former mass corresponds to the lightest right-handed neutrino 
mass induced by the loop effect, although the latter one to the SM gauge singlet fermion. We 
show each of the DM annihilation processes and compare to the observation of relic abundance, 
together with the constraints of Lepton Flavor Violation (LFV) and active neutrino masses. 
Moreover we show that our model has some parameter region allowed by the direct detection 
result reported by XENON100, and it is possible to verify the model by the future XENON 
experiment. 
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1 Introduction 



It has been verified that neutrinos have tiny masses by neutrino oscillation experiments [IH7]- Unfor- 
tunately the finite neutrino masses are not explained in the framework of the Standard Model (SM). A 
lot of models have been proposed to extend this point [8HTTj . On the other hand the existence of non- 
baryonic Dark Matter (DM), which dominates about 23% of the Universe from the CMB observation 
by WMAP [12], is also shown by the cosmological observations in our Universe [13j[T3] 0. More- 
over direct detection experiments of DM are performed around the world such as XENON100 |16j . 
CRESSTII [IT], CoGeNT [IS], DAMA [IS] and TEXONO |2Q]. Especially XENON100 experiment 
gives the most severe limit for elastic scattering cross section between DM and nucleon [16]. This 
implies that DM in the Universe interacts very weakly with quarks. It would be that DM has no 
interaction with quarks. DM is required in the Universe, however a candidate particle of DM is also 
not included in the SM. Although the property of neutrino in the SM is similar with that of DM, 
neutrinos are too light to be DM candidate. Thus in order to improve this problem, it is necessary to 
add new particles as DM candidate in the SM. Therefore, these current experiments about neutrinos 
and DM suggest serious verifications that the SM should be modified in order to accommodate the 
existence of DM as well as non- vanishing neutrino masses. 

Radiative seesaw models are known as attractive framewoks for new physics at TeV scale that 
can provide an elegant solution to explain these two matters of grave concern simultaneously [2TH21] . 
This kind of model correlates the finite neutrino masses with the existence of DM since neutrino 
masses are generated by radiative effect and DM runs inside the loop. In particular, the radiative 
seesaw model proposed by Ernest Ma [21] is one of the simplest models. Subsequently there are a lot 
of recent works in terms of the model [251427] and the extended models [28HH] . The other models of 
radiative neutrino mass are studied in Refs. jl2] - [5"T] . 

In this paper, we propose a new model of two-loop induced neutrino masses with local B — L 
symmetry. Due to the two remnant Abelian symmetries (Z2 and Tjq) even after the B — L and 
electroweak symmetry breaking, our model has multi-component DMs. Two or three particles of 
them can be DMs simultaneously depending on the mass hierarchy. Moreover since one of DMs 
obtains the mass at one loop, we expect it to be rather light with mass of (9(10) GeV. We check 
whether they can satisfy the correct relic density of DM observed by WMAP, and also the upper 
bound of elastic cross section with nucleon by XENON100. 

This paper is organized as follows. In Section 2, we show our model and discuss the Higgs sector 

1 Very recently, new result of the DM relic density was given by Planck measurement as Q,c.h 2 = 0.1196±0.0031 [15j . 
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including the Higgs potential, stationary condition, S-T parameters and neutrino mass in lepton 
sector. In Section 3, we analyze DM phenomenologies. We summarize and conclude in Section 4. 

2 The Two-loop Radiative Seesaw Model 
2.1 Model setup 
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Table 1: The particle contents and the charges for fermions. 
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Table 2: The particle contents and the charges for bosons. 

We propose a two-loop radiative seesaw model with U(1)b-l which is an extended model of 
the seesaw model of Ma [21] . The particle contents are shown in Tabs. Q] and [2j We add three 
right-handed neutrinos N c , three SM gauge singlet fermions S and S, a SU(2)l doublet scalar t] 
and B — L charged scalars x an d £ to the SM content, where t] and \ are assumed not to have 
vacuum expectation value (VEV). The B — L charged scalar S is the source of the spontaneous 
B — L breaking by its VEV of (£) = v'/V2 ~ O(10) TeV. The Z 2 parity is also imposed so as to 
stabilize DM candidates. The right-handed neutrinos N c do not have masses at tree level. As a 
result, the neutrino mass is obtained not through the one-loop level (just like Ma-model [21]) but 
through the two-loop level. 
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The renormalizable Lagrangians for Yukawa sector and Higgs potential are given by 

-Yukawa = ye&e c L + y v rfN c L + y N N°xS + ys^SS + y§TtSS + h.c, (2.1) 
-£ H igg S = ml&® + mlrfrj + m^E^E + m\x^X + m s[x 2 E + h.c.] 

+A 1 ($ t $) 2 + X 2 {rfr)) 2 + A 3 ($ t $)(r ? t r ? ) + X 4 (&rj){rf$) + \ b [{&r)) 2 + h.c] 
+A 6 (£t£) 2 + A 7 (StS)($t$) + A 8 (E+E)(^) + Ux^x) 2 

+Aio(x t x)($ t $) + An(x f x)(^) + A 12 |yS| 2 , (2.2) 

where A5 has been chosen real without any loss of generality. The couplings Ai, A2, A6 and A9 have to 
be positive to stabilize the Higgs potential. Inserting the tadpole conditions; m 2 = — A^ 2 — \ 7 v' 2 /2 
and m\ = —Xqv' 2 — Xyv 2 /2, the resulting mass matrix of the neutral component of $ and E defined 

as 

o _ v + 0°(x) _ v' + a(x) 

is given by 

2 ^ / 2Aif 2 Xfuv' \ f cos a sin a \ / ml \ / cos a —sin a \ 
V A7W 2A6f /2 y y — sin a cos a J y m 2 ^ y y sin a cos a y 

where /i implies SM-like Higgs and if is an additional Higgs mass eigenstate. The mixing angle a is 
given by 

tan 2c* = X 2 7m i a - ( 2 - 5 ) 
The Higgs bosons <p° and a are rewritten in terms of the mass eigenstates h and if as 

c/>° = /i cos a + if sin a, 

a = —h sin a + if cos a. (2-6) 



The other scalar masses are found as 



l.o 1 



m\=m 2 {jf) = m l + lX z v 2 + ^XW\ (2.7) 
m 2 R =m 2 (Rer ] °) = m\ + ^A 8 ^ + ^(A 3 + A 4 + 2A 5 > 2 , (2.8) 
m 2 7 =m 2 (lmr]°) = m 2 2 + lx 8 v' 2 + ^{X 3 + X 4 - 2X 5 )v 2 , (2.9) 



2 2 

L = m 2 + ^A 10 ^ 2 + ^X 12 v' 2 + V2m 5 v', (2.10) 
1 1 

m 2 = m 2 A + -X w v 2 + -X 12 v' 2 - V2m 5 v'. (2.11) 
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The tadpole conditions for r\ and x, which are given by 



av 



VEV 



0. 



9V 



VEV 



0, < 



a 2 v 

drj' 2 



and 



< 



d 2 V 
&x 2 



VEV 



tell us that 



VEV 
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0<^2 + y(^3 + A 4 + 2A 5 ) + yA 8 , < m\ + V —\ w + \/2m h v' + yA 12 , (2.12) 

in order to satisfy the condition v v = and w x = at tree level, respectively. In order to avoid that 
(<&) = (S) = be a local minimum, we require the following condition: 
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To achieve the global minimum at (rj) = (%) = 0, we find the following condition 



(2.13) 



< An - -a/A 2 A 9 . 



Finally, if the following conditions 



< A 3 + - VA1A2, 0<A 7 + -\A 1A6, < A10 + -a/ AiAc 



0<A 8 + - v / AA, < An + -VA 2 A 9 , < Ai 2 + -v/A^Ag", 



(2.14) 



(2.15) 



are satisfied, the Higgs potential Eq. fl2.2p is bounded from below. 



2.2 S and T parameters 

It is worth mentioning about the new contributions to the S and T parameters due to the new scalar 
boson rj, which are given in Refs. j52j|53] as 
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b 2 \" b 2, 

where a cm = 1/137 is the fine structure constant. The experimental deviations from the SM predic- 
tions, under m hsM = 126 GeV, are given by [5 



0.03 ±0.10, T T 



0.05 ±0.12, 



(2.19) 
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Figure 1: Neutrino mass generation via two- loop radiative seesaw. 



When the masses are 1 < mi/rri2 < 3, the function F(m 1 ,m 2 ) is approximated to 

2 2 
F(mi,m 2 ) ~ - (mi - m 2 ) , 

as in Ref. [52]. From Eq (I2.19p . we get the following constraint for r\ masses, 



{m v - m VR ) (m v -m m )< 133 GeV. 



(2.20) 



(2.21) 



2.3 Neutrino mass matrix and LFV processes 

The active neutrino mass matrix at the two-loop level as depicted in Fig{T] is given by 

m s r 1 . r 1 -* , i 
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(2.23) 
(2.24) 



x(l — x) 

and 777,5 is the mass of S, abbreviating generation index of S. The neutrino mass scale should be 
roughly m v ~ 1CT 1 eV. This implies that Yukawa coupling is around y 2 N yl ~ 1CT 8 when 7775 ~ 1 TeV 
and the integration by x and y to be of 0(1) are assumed. 

The Branching Ratio (Br) of charged Lepton Flavor Violation (LFV) processes £ a — > £37 (a,j3 = 
e, fx, t) is given by 



Br {£ a ->£fn) 
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where right-handed neutrino mass is neglected here. The latest limit for [i — >■ cy is given by MEG 
experiment [SS] as 

Br(/i e7) < 5.7 x 1(T 13 . (2.26) 

For sum of active neutrino masses, the limit of Yl m v < 0.933 eV is imposed from the cosmological 
observation [15]. In the next section, we take into account these constraints of S-T parameters, LFV 
and the neutrino mass in the discussion of DM. 



3 Dark Matters 

We discuss the DM properties in this section. The Z 2 parity imposed to the model stabilizes DM. 
In addition to the Z 2 parity, we have a remnant Z 6 symmetry after B — L symmetry breaking which 
stabilizes particles charged under the Zq symmetry as well. As a result, our model has two or three 
DM candidates simultaneously, and the number of DMs depends on mass hierarchy of DM candidates 
included in the model. In general we have five DM candidates which are N c , S, S as fermionic DMs 
and Xr{i)i Vr(i) as bosonic DMs (For bosonic part, the DM proparty of the imaginary part is almost 
same as that of the real one). Of these, three particles can be DMs when decay of a charged particle 
under the Z 2 and/or Zg symmetry is kinematically forbidden, otherwise we have two DMs. The mass 
of N c is expected to be somewhat light (~ (9(10) GeV) because its mass is generated at the one-loop 
level. The other particles have typically the mass of B — L symmetry breaking scale. Therefore it 
is natural to choose N c as the lightest DM. In this case, rjn and rjj cannot be DM since they have 
the same charge with N c under Z 2 and Zg symmetry. The remaining DM candidates are S, S, xr 
and xi- The interactions of S and S are almost same except for Yukawa interaction y^N c xS. This 
Yukawa interaction leads DM-exchanging scattering like SS — > N C N C . The DM candidates xr an d 
Xi have also DM-exchanging scattering Xr(i)Xr(i) ~~ >" N C N C via Yukawa interaction yNN c xS as same 
as the case of S. Thus as the simplest case, we consider two-component DMs (N c and S) in the 
following; since there is no exchange scattering between W « a,d Sat the tree levelg 

The mass matrix of three right-handed neutrinos N c is radiatively induced and the expression is 
found as 



(m NC ) - y (yN)ik(yN) j km Sk r ™ 2 XR ^ / m\ R \ m 2 Xi ^ / m? Xl \ 



(3.1) 



2 Regarding the pair of (N C ,S) or (N c ,xr(i)) as DMs. it would be the simplest if Yukawa coupling ys is large 
enough because most of S or XR{i) annihilates into N C N C . In this case, one can consider that only 7V C is substantial 
DM. 



6 




h 



f s 



w-/z s 



Figure 2: The annihilation channels of fermionic DM N (upper row) and S (lower row) 



We need multi-pair of S and S to generate three non-zero right-handed neutrino masses, otherwise 
only one non-zero mass is generated. An interesting feature of the model is the connection between 
the light neutrino masses Eq. (I2.22p and the right-handed neutrino masses Eq. (13.11) . The right- 
handed neutrino masses tend to be small corresponding to the tiny light neutrino masses. The mass 
matrix can be diagonalized by a unitary matrix and the lightest particle, which is called as simply 
N below, is DM. We choose the diagonal basis of right-handed neutrinos. The annihilation channels 
of the right-handed neutrino DM into the SM particles are shown in the upper part of Fig. [21 where 
the annihilation channel via B — L gauge boson is omitted since the contribution is small enough 
due to exchange of heavy B — L gauge boson. We obtain the annihilation cross sections as 

2 
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(3.4) 



where G (x, y) = xy (1 — x — y + 2xy), a v = m 2 N / (m 2 N + rrvt) and a a = m 2 N / (m 2 N + m 2 a ), a = R, I. 
In the last expression, the symmetric factor 1/2 should be multiplied when the flavor of the final 
state vlVl is the same. The chiral suppression is not effective for the pair of (anti-)neutrino final 
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state channels if the mass splitting between r)R and rji is not negligible. 

We analyze the relic abundance of the lightest right-handed neutrino N with constraints from 
S-T parameters, the neutrino mass and LFV Eq. (12.261) . We sweep the parameters in the following 
range: 



The mass of r\ must satisfy the constraint of S-T parameters Eq. (12.211) . and the typical mass scale of 
S, x an d tj is of 0(1) TeV since we assume the B — L breaking scale is several TeV. Yukawa couplings 
(y v ) Tl and yN are fixed to {y u ) Tl = Un = 1-0 in order to have a proper annihilation cross section, 
being consistent with /i — > cy. One might think that Yukawa coupling {y v ) T i = 1 conflicts with 
the scale of active neutrino masses. However it can be kept to the appropriate scale by considering 
structure of y^ appearing in Eq. (I2.22p . As concerning above factors, the right-handed neutrino 
masses shown in Eq. (13.1 p induced by the one-loop effect become (9(10) GeV, otherwise the right- 
handed neutrino masses become much lighter and the annihilation cross section will be too small 
because it is proportional to a Vj Rj. 

The result in the case of N DM is shown in Fig. [31 and each point satisfies the constraint from 
S-T parameters, appropriate light neutrino mass scale ~ 10" 1 eV, the LFV constraint Br (/i — > ej) < 
5.7 x 1CT 13 and thermal DM relic density which corresponds to <™ r ei > 3 x 10" 26 cm 3 /s. The reason 
of the inequality in the cross section is that we have two DM candidates. One should note that 
total relic density is supplied by N and the other DM S. Therefore even if the annihilation cross 
section is larger than 3 x 10~ 26 cm 3 /s, the lack of the amount of DM is compensated by the amount 
of S. One can see from the upper left panel that most of parameter space can satisfy the current 
upper bound of the branching ratio for fi — > &y. This is because Yukawa couplings relating with 
jj, — > e7 are the order of 10~ 3 ~ 10 -2 as shown in the upper right panel. The right-handed neutrino 
mass should be of (9(10) GeV in this case. The magnitude of Yukawa couplings is also important to 
obtain the proper light neutrino mass scale (See Eq. (12.221) ). In addition to the Yukawa coupling, a 
small mass splitting is somewhat required to obtain correct light neutrino masses as the left lower 
panel shows. As a result, our scenario of DMs is still p-wave dominant since Eq. (13.41) is suppressed. 
Due to the small 7] mass splitting, we can get tiny neutrino masses because of a cancellation among 
the integrands of Eq. (12.221) . The masses of rjR and t]i are restricted to be roughly 10 2 GeV. Thus 
charged rj mass is around 200 GeV from the constraint of S-T parameters. On the other hand as 
shown in the right-lower panel, mass hierarchy between \r an d Xi is necessary to get the scale of 



10 2 GeV <m s < 10 4 GeV, 10~ 3 < (y„) el 
10 2 GeV < m v < 10 3 GeV, 10 2 GeV < m 
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n < 10 4 GeV. 
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Figure 3: The parameter spaces of satisfying LFV, the light neutrino mass scale, the thermal DM 
relic density. 



the right-handed neutrino masses which is connected with the size of the annihilation cross section. 

Next, we move on discussion of S DM. This DM does not have any interactions with N at tree 
level. Hence we can consider these two DMs separately. The annihilation cross sections of S shown 
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in the lower part of Fig. [2] are given as 
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where mg is the mass of S and the color factor Cf is 1 for leptons and 3 for quarks. The parameter 



is defined as 
Higgs H, 



,/ (2m| — m^) and -D is the propagator of the SM-like Higgs h and an extra 
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(3.11) 



The if if final state process can be obtained by replacing rrih — > ran and sin a — > cos a in Eq. (I3.10p . 
We have m<j, ?/<j, sin a, m# as parameters and sweep in the following range: 



200 GeV < m S(H) < 5000 GeV, 10" 3 < y s < 1, 10" 3 < sina < 1, 



(3.12) 



satisfing <7t> > 3 x 10~ 26 cm 3 /s. The result is shown in Fig.HJ The left panel shows that the mg < mjj 
roughly holds when rrig is larger than around 1 TeV. Moreover the couplings have to be of order one 
to get correct relic density of S (the right panel). It causes the elastic cross section with nuclei to be 
larger than the upper bound by XENON100 as we will discuss below. 

Let us move on to the discussion of direct detection of DMs. Although our DM consists of two 
components iV and S, N does not have any interactions with quarks at tree level since it is a right- 
handed neutrino! 3 ] The other DM S interacts with quarks via Higgs exchange. Thus it is possible 
to explore the DM in direct detection experiments like XENON100 [16]. The Spin Independent (SI) 
elastic cross section ogi with proton p is given by 
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(3.13) 



3 The DM N can interact with quarks at loop level through electromagnetic interactions [25l[56] and it would be 
large as same as being detected by the XENON experiment. 
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Figure 4: The allowed parameter spaces for S DM. 



where fig = (rn^ 1 + is the DM-proton reduced mass. The parameters which imply the 

contribution of each quark to proton mass are calculated by the lattice simulation [5*71 11)8] as 

ft = 0.023, f A = 0.032, f v s = 0.020, (3.14) 

for the light quarks and fg = 2/27 ^1 — X^<3 fq^ f° r the heavy quarks Q where q < 3 implies the 
summation of the light quarks. The recent another calculation is performed in Ref. [59|. Fig. [5] 
is the comparison with XENON100 upper bound (16] with the same parameter obtained from the 
analysis of the relic density. In the case that S DM to be dominant, most of parameter region 
allowed by the relic density is excluded by the XENON100 upper bound due to the large Yukawa 
coupling y s . Despite of such a strong constraint, some allowed parameter region certainly exists. 
These parameters imply that Yukawa coupling y s (with large mixing of a) is rather small and the 
mass of m§ is close to a resonance for the annihilation cross section in Eq. fl3.7p - fl3.9p . Needless to 
say, we can easily relax such a situation because of multi-component DM scenario. Since we have 
two DMs in the model, the fraction parameter of relic density £, which stands for the fraction of 
relic abundance of S to the total abundance, makes the XENON100 limit looser, and wide allowed 
parameters appear. 



11 



XENON100 (2012) (£=1.0) 
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10 ' ' ' ' ' — J — ' ' ' 
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Figure 5: The comparison of the elastic cross section of S with nucleon where the XENON100 upper 
bound is also drawn together where the parameter £ stands for the fraction of relic abundance of S 
to the total abundance. 

4 Conclusions 

We have constructed a two-loop radiative seesaw model with local B — L symmetry at the TeV scale 
that provides neutrino masses. We have also studied the multi-component DM properties, in which 
we have two fermionic DMs with different mass scale; (9(10) GeV for N and O(10 2 ~ 3 ) GeV for S. 
Although N is right-handed neutrino, its mass is generated by the one-loop effect. We have shown 
that the allowed mass regions of the particles t)r(i), Xr{i)i an extra Higgs H, S and Yukawa couplings 
constrained by S-T parameters, Br(/^ — > e j)< 5.7xl0~ 13 , m v ~ 0.1 eV, and annihilation cross 
section of DMs, in which for example we found the mass of t)r and rji should be softly degenerate: 
\m r)R — m Vl \ — lO - ^ 1-2 ' GeV. The upper bound of is around 30 GeV, which is naturally satisfied 
because of the loop-induced mechanism of m^. Too light leads too much relic density of N 
because of the small annihilation cross section. Verification of our model is also possible by direct 
detection of S DM through the interaction with Higgses. In particular, the region of the large mixing 
sin a will be testable by the exposure of the XENON experiment. 
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